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ABSTRACT 

We study the stellar angular momentum of thousands of galaxies in the Illustris cosmological simulation, 
which captures gravitational and gas dynamics within galaxies, as well as feedback from stars and black holes. 
We find that the angular momentum of the simulated galaxies matches observations well, and in particular 
two distinct relations are found for late-type versus early-type galaxies. The relation for late-type galaxies 
corresponds to the value expected from full conservation of the specific angular momentum generated by 
cosmological tidal torques. The relation for early-type galaxies corresponds to retention of only ~ 30% of 
that, but we find that those early-type galaxies with low angular momentum at z = 0 nevertheless reside at 
high redshift on the late-type relation. Some of them abruptly lose angular momentum during major mergers. 
To gain further insight, we explore the scaling relations in simulations where the galaxy formation physics is 
modified with respect to the fiducial model. We find that galactic winds with high mass-loading factors are 
essential for obtaining the high angular momentum relation typical for late-type galaxies, while AGN feedback 
largely operates in the opposite direction. Hence, feedback controls the stellar angular momentum of galaxies, 
and appears to be instrumental for establishing the Hubble sequence. 

Subject headings: galaxies: formation — structure — fundamental parameters — kinematics and dynamics — 
methods: numerical — hydrodynamics 


1. INTRODUCTION 

The angular momentum content of a galaxy, which corre¬ 
lates with galaxy mass and type, is one of its primary proper¬ 
ties. The cosmological structure formation framework pro¬ 
vides th e basi c tenet of tidal torque theory (iPeeblesl 119691 
iDoroshkevichl 1 197(1) . Within the ACDM paradigm, the an¬ 
gular momentum J of dark matter halos, expressed using the 
‘spin parameter’ A = / \f~E/GM 5 ^ , with E representing en¬ 
ergy and M mass, is robustly measured in N-body simu lations, 
having a typical value of A ~ 0.034 (e.g. iBett et all 120071) . 
However, deriving the angular momenta of galaxies requires 
further assumptions. In the ‘standard’ theory of disk galaxy 
formation, baryons retain their specific angular mom entum 
/ = J/M as_they collapse to halo centers dFall & Efstathioul 
119801 iMo et al.lll998l) . There are, however, some reasons to 
question the generality of a simple relation between the angu- 
lar momenta of galaxies and their dark halo s (Ivan den Boschl 
l200HlVitvitska et al.ll200l IBett et al.ll20Toh . 

With these issues in mind, we define an ‘angular momen¬ 
tum retention factor’ as the ratio between the specific an- 
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gular momenta of a galaxy’s stars and dark matter, r/ ; - = 
j* /./dm - The retention factor should be considered in light 
of two factors, (i) The baryons ending up in the galaxy may 
not uniformly sample the initial angular momentum distri¬ 
bution of all the baryons . Thi s can arise from th e accre¬ 
tion process (iKassin et alJl2012t IStewart et all 120131). which 
is complex and far from s pherically-symmetric dKeres et all 
l200l INelsonetaDl20l5ah. as well as from gas expulsio n 
by galactic outflows (iBinnev et al.M2Q01b IBrook et all 1201 11) . 
(ii) The angular momentum may be modified with respect to 
its initial value. Losses can occur on the way to the galaxy 
( GanetalJ 120101: iDanovich et all 1 20141), d uring mergers 
( Hernquist & Mihosll 1995 1: Jesseitet^ 20091). and inside th e 

dDekel et all I 


galaxy itself (iDekel et all|2009h iHonkins & Q uataertj|2011|). 
Gains can be induced by galactic fountains (IBrook et al.l2012l : 
lUbleretal.ll20Ti) . 

The most general tool to study these complexities of an¬ 
gular momentum buildup towards the overall rjj is cosmo¬ 
logical hydrodynamical simulations. Those simulations self- 
consistently follow the baryonic non-linear dynamics from 
cosmological initial conditions, and allow modeling a va¬ 
riety of strongly-coupled galactic processes. Early cosmo¬ 
logical hydrodynamical simulation s su ffered from a ‘catas¬ 
trophic’ angular momentum loss dNavarro et al .111 993 ). Later 
work sh owed reduced losses a s a result of improved res¬ 
olution dGovernato et akl [20071: iKa ufmann et al.l 20071) . im¬ 
proved numerical techniques (ISiiacki et all 20121), and in¬ 
clusion of stronger f eedback dSommer-Larsen et all Il999t 
iGovernato et al.ll2007 :j. Recent simulations started exploring 
galact ic angular momentum distinguishing different galaxy 
types dNaab et al.ll2014l : lFiacconi et al.ll2014l) . Only the most 
recent generation of cosmological hydrodynamical simula¬ 
tions reproduces a variety of galaxy types in large numbers, 
allowing a more detailed andjnore favorable comparison to 
observations (e.g. lTeklu et al.ll2015l) . 
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An early account of observed galaxy angular momenta 
dFalll 119831) . confirmed by more sophisticated later anal¬ 
ysis, identified that angular momentum depend s bo th on 
galaxy mass and mor phological type dRomanowskv & Falll 

2012: lFall & Romanowskv1l2013l[Qbreschkow & Glazebroold 

201 4). First, higher-mass galaxies have a higher specific an¬ 
gular momentum. Second, early-type galaxies contain about 
5 times less angular momentum than late-type galaxies of the 
same stellar mass. Moreover, ij j appears to be roughly mass- 
independent for each of the two main galaxy types. Disk 
galaxies have r/j ~ 0.8, while early-type galaxies on average 
lie on a parallel scaling relation characterized by r/ ; ~ 0.1-0.2 
dFall & Romanowskvi2013l) . 

In this letter, we study the angular momentum of thou¬ 
sands of galaxies in a set of cosmological hydrodynami- 
cal s imulations. Our fiducial r un is the Illustris simula¬ 
tion (iVogelsberger et al~l I2014allbl : iGenel et ahl 120141) . which 
uses a physical model that results in qu ite a g ood mat ch 
to a set of observed sc aling relations dVbgelsberg er et alJ 
1201 3b iTbrrev et all 120141) . In particular, a qualitative agree¬ 
ment with observed mass- and redshift-dependent trends was 


found for galaxy mori 

Dhologies and colors (Vogelsberger et al. 

2014a,ti Genel et al. 2014 Torrey et al. 20151 Snyder et al. 

2015|). This paper is 

structured as follows. In Section [2] we 


briefly describe the simulations and analysis methods. In Sec- 
tion[3]we study the angular momentum content of galaxies in 
our simulations, and in Section [4] we summarise our results 
and conclude. 


2. METHODS 

The primary si mulation used in t his pa pe r is the Il¬ 
lustri s simulation C Vogelsber ger et al.l l2014allbl : IGenel et al J 
120141) . which evolves a ACDM cosmological volume of 
(106.5 Mpc) 3 to z = 0. Gravity and hydrodynamics, as well 
as additional galaxy formation processes, are_modeled us¬ 
ing the AREPO TreePM-moving-mesh code dSpringelll2010l) . 
With « 2 x 1820 3 resolution elements, the typical baryonic 
particle mass is ps 1.26 x 10 6 Mq, and the gravitational soft¬ 
ening is 1.4kpc comoving (capped at 0.7 kpc physical at 
z = 1 for baryonic particles). The code implements phe¬ 
nomenological, parameterized models for gas cooling, star 
formation and the resulting galactic winds, metal produc¬ 
tion and mass return driven by stellar evolution, and black 
hole f ormation, accretion , and fe edback, all described in de¬ 
tail in IVogelsberger et~ahl (12013!) . We calibrated the fiducial 
model such that the history of star formation density and 
the z = 0 relation between stellar mass and halo mass are in 
reasonable agreement with observations (iVogelsberger et all 
120131) . Notable tensions do remain, however, with over¬ 
production of galaxies both above and below the ‘knee’ 
of the mass function. Following this calibration, we have 
made extensive comparisons to other observables, finding 
broad agreement in many p arame ters the_model was not 
tuned for fe.g.lBir d et al.l2014 :lRodriguez-G omez et all2015t 

ISales et akll2015USnvder et alll2015HWellons et al.ll2015lf 

In addition, a smaller volume of (35.5Mpc) 3 has been 
evolved with variations around the fiducial model. In some 
runs certain model ingredients were turned off, and in oth¬ 
ers the st rength of certain feedback processes were d ialed up 
or down (IVogelsberger et al1l201 3t iTorrev et akll2Q 141) . These 
simulations have a lower resolution, using a gravitational soft¬ 
ening of 2.8 kpc comoving (capped at 1.4kpc physical at 
z = 1 for baryonic particles), and baryonic particle masses of 


~1.7x 1O 7 M 0 . 

G alaxi es are identified using the SUBFIND algorithm 
(ISnringel et al.l |200H) . We include in the angular momen¬ 
tum calculation all stellar particles that belong to the subhalo, 
i.e. those that are bound to its primary density peak but do 
not belong to its satellites. Herein we apply no restriction on 
radius, since a significant fraction of the angul ar mom entum 
usually lies beyond the half-mass radius dRomanowskv & Falll 
120121) . However, we verified that restricting the measure¬ 
ment to a radius five or ten times the stellar half-mass ra¬ 
dius, roughly correspondin g to the observational estimates 
dRomanowskv & Fall 12012ft . would not change our results 
substantially (but a stronger re strictio n on radius would be¬ 
come significant, see also lTeklu et al.ll20l5l) . 

For the calculation, each galaxy is centered on the particle 
with the lowest gravitational potential of the subhalo. It is 
important not to use the center of mass position, since it is 
sensitive to structure at large radii, and may not represent the 
rotation center of the galaxy. However, the velocity of the cal¬ 
culation frame is set to move with the stellar center of mass, 
which is, importantly, not sensitive to the velocity of any in¬ 
dividual particle. 

3. RESULTS 

Each panel of Fig. [j] shows the distributions (gray), and 
their means (black), of stellar specific angular momentum 
j * as a function of stellar mass M*, for all Illustris galax¬ 
ies with M* > 10 9 Mq at z = 0. These distributions are uni- 
modal. However, the means for two ‘extreme populations’, 
selected based on a third, independent quantity - flatness, 
concentration, or star-formation rate (left to right, respec¬ 
tively) - show a clear separation (cyan, purple). These are 
the 10% tails on each end of the distributions of those pa¬ 
rameters. Selecting wider tails of 30% at each end results 
in factor ps 2 reduction of the logarithmic separation, with 
variation between panels. The left-hand panel selects galax¬ 
ies based on axial ratio. To this end, we center the galaxy 
on the most bound particle and align the z axis with the an¬ 
gular momentum vector. Then we measure the second mo¬ 
ment of the mass distribution in each cartesian direction, 
Mj = (J2 m krki) 1 / 2 , where k enumerates over stel- 

k ’ k 

lar particles, and i £ (x,y,z). We define ‘flat’ galaxies as those 
having M z /(M x M y ) l/2 < 0.55, and ‘round’ galaxies to have 
M z /(M x My ) 1 / 2 > 0.9. In the middle panel, we separate galax¬ 
ies based on their morphology measured in the unattenuated 
; band using the Gini-AUo structural parameters dLotz et al.l 
12004 ). Specifically, we define the ‘ F(G,M 2 q ) bulge statis¬ 
tic’ ( Snyder et al.l 1201 51) and compute its average over four 
random viewing angles for each galaxy. We define ‘diffuse’ 
galaxies as those with F < —1 (small bulges) and ‘concen¬ 
trated’ galaxies as those with F > 0.1 (significant bulges). 
In the right-hand panel, galaxies are separated by specific 
star-formation rate, with ‘star-forming’ galaxies satisfying 
sSFR > 0.2 Gyr -1 , and ‘quiescent’ ones sSFR < 0.01 Gyr -1 . 
Finally, the right -hand panel also shows two visually-selected 
sets of galaxies (IVogelsberger et akll2014bl) . 42 ‘blue’ galax¬ 
ies and 42 ‘red’ ones, a selection that correlates strongly with 
visual morphology. 

In addition to the simulated resul ts, Fig .[pshows fits to the 
observed correlations (dashed; iFall & Romanowskvl 120131) . 
for late-type (blue) and early-type (red) galaxies. While we do 
not have a visual classification for each of the 29,203 Illustris 
galaxies on this figure, the three cuts shown in the different 
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Figure 1 . The ; = 0 relation between stellar specific angular momentum and stellar mass. In all panels, dashed curves show the IFall & Romanowskvl 12013 !;) 
observed relations, in the range where data exist, for disks (blue) and spheroids (red), which both have a slope of R3 2/3. The black curves show the logarithmic 
mean relation for all Illustris galaxies, whose full distribution appears in gray. Each panel shows two additional mean relations for ‘extreme’ galaxies, selected 
to have high (cyan) or low (purple) values of an independent quantity. Plotting median values, rather than mean, results in essentially indistinguishable curves. 
Left: Separation by axial ratio. Middle: Separation by morp hology. Right: Separ ation by specific star-formation rate, and in addition the blue (red) squares show 
the visually-selected late- (early-)type Illustris galaxies from IVogelsberger et all <2014b :i. 


panels represent reasonable proxies of late- versus early-type 
galaxies. Each of these cuts results in a somewhat different 
‘angular momentum sequence’, but as the late-type (cyan) and 
early-type (purple) relations separate by factors of a few in 
most cases, overall they clearly demonstrate the correlation 
between galaxy type and specific angular momentum in Illus¬ 
tris. The late-type relations are somewhat shallower than the 
observed one, and generally have lower normalizations, while 
the slope and normalization of the early-type relations match 
observations well. In two of the three cases (left and right pan¬ 
els), the relation for the full galaxy population (black) is closer 
to the late-type relation at low masses, and to the early-type 
relation at high masses, representing the mass-dependent pro¬ 
portions between the two galaxy types. Keeping in mind that 
the uncertainties in the observed quantities, both angular mo¬ 
mentum and stellar mass, are about ±50%, we conclude that 
galaxies at z = 0 in Illustris are consistent with having angular 
momentum content in good agreement with observations. 

In Fig. [2]we show the same data, except we replace the stel¬ 
lar mass on the x-axis with the dark matter halo mass of each 
galaxy. Here we only include central galaxies, and note that 
Fig- HI itself would only change negligibly had we excluded 
satellite galaxies there as well. This allows us to compare 
the angular momentum content of the galaxies to the angular 
momentum they are expected to have under the assumption 
of angular momentum conservation during galaxy formation, 
and a halo spin parameter of A = 0.034 (brown). This compari¬ 
son shows that ‘late-type’ galaxies (as defined by the different 
cuts and represented by the cyan curves) have an angular mo¬ 
mentum content very close to that with no angular momentum 
loss, i.e. rjj ~ 100%. The ‘early-type’ galaxies form a paral¬ 
lel sequence corresponding to only r/ ; ~ 30%. Two puzzles 
immediately emerge from this result, (i) Have the baryons 


forming the late-type galaxies indeed avoided losing the an¬ 
gular momentum they acquired by tidal torques when they 
were still in intergalactic space, or do other mechanisms con¬ 
tribute to shaping their sequence? (ii) What is the physical 
cause of the mass-independent mean r/j value of ~ 30% (or 
loss factor of ~ 70%) of early-type galaxies? 

To probe the first of these questions, we explore in Fig. [3] 
the consequences of modifying the feed back sub-grid mod¬ 
els included in the simulation. Fig. |3(a)| presents the specific 
angular momentum versus halo mass of central galaxies in 
a series of cosmological simulations of a (35.5Mpc) 3 vol¬ 
ume with increasing completeness of the physics models. We 
first observe that galaxies in a simulation without any signif¬ 
icant feedback (green) have typically retained only ~ 30% of 
their angular momentum with respect to the tidal torque value 
of A = 0.034 (see also iTorrev et all[2012l) . and hence have 
a typical angular momentum content close to that observed 
for early-type galaxies. Adding mass return from evolving 
stellar populations and metal-line cooling (beige) does not 
change the angular momentum content. Introducing galac¬ 
tic winds (orange), however, increases the angular momen¬ 
tum dramatically across all ma ss scales. The effect of galactic 
winds is seen also in Fig. |3(b)[ where the thick (thin) solid or¬ 
ange curves show simulations with increased (reduced) mass¬ 
loading factors for the galactic winds (by factor 2), indicating 
that stronger winds ten d to produce galaxies with higher angu¬ 
lar m om enta (see also [Sc annaniec o et alJ 120081 IZavala et alJ 
120081 but IS ales et alJl2Q10l for a contradictory result). 

Importantly, in the presence of galactic winds, the angu¬ 
lar momentum increases to a level approximately consistent 
with full retention of angular momentum. However, we spec¬ 
ulate that this may not indicate strict conservation of angu¬ 
lar momentum, but rather a combination of loss on one hand 
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Figure 2. The relation between stellar specific angular momentum and dark matter halo mass (enclosing a density 200 times the critical density). The various 
Illustris data are the same as in Fig. [7] except for the omission of satellite galaxies. The specific angular momenta of halos with a spin parameter of A = 0.034 
is displayed for reference (brown). The results from the two lower-resolution versions of the Illustris simulation (not shown) deviate from these results by up to 
0.1 dex for late-type galaxies, while for early-type galaxies they are very well converged. 



Figure 3. The z = 0 relation between stellar specific angular momentum and halo mass in simulations with different physics models. The specific angular 
momentum of halos with a spin parameter A = 0.034 is displayed as well (brown). Left: Increasing the completeness of the model, from no feedback up to the 
full model. Right: Modifying the total power of galactic winds (solid orange) or radio-mode feedback (dashed cyan) by a factor of 2. We find that galactic winds 
increase the angular momentum of galaxies, while AGN radio-mode feedback reduces it. 
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and gain on the other by the presence of a strong ‘halo foun¬ 
tain’. In our simulations, gas that is ejected into the galac¬ 
tic winds typically comes back to the galaxy after spending 
some time in the halo, and this happens many times over 
dMarinacci et al.l2014tiNelson et al.l2015bl) . The gas can then 
acquire angular momentum in the halos, and also induce halo 
gas with high angular momentum to cool down to the galaxy 
dBrook et al.ll2012HUbler et~aT]l2014l) . We leave a quantitative 
analysis of this complex exchange process to future work. 

Fig. |3] further demonstrates that quasar-mode AGN feed¬ 
back (left panel, cyan) hardly modifies galactic angular mo¬ 
mentum. However, adding, or enhancing, the radio-mode 
AGN feedback (where thermal bubbles are injected into halo 
atmospheres), has the opposite effect to winds, i.e. it reduces 
the angular momentum content of galaxies. This can be un¬ 
derstood in terms of the effectiveness of this type of feed¬ 
back at late cosmic times. Gas accreted onto halos at low red- 
shift has a higher specific angular momentum, but its accre¬ 
tion onto the galaxies is suppressed by a stronger radio-mode 
feedback, leaving the galaxies with preferentially more early- 
accreted baryons, which have lower specific angular momen¬ 
tum. 

Next we probe the second puzzle presented earlier, namely 
why early-type galaxies show a mean r]j that is independent of 
mass, and thereby form a parallel relation_to that_of late-type 
galaxies. We use merger trees dRodriguez-Gomez et al.l2015h 
to present in Fig. |4l the evolution tracks of vi sually-selected 
early-type galaxies (1 Vogelsbersier et al1l2014bl) in specific an¬ 
gular momentum and stellar mass, from z = 3 to z = 0. We 
find that early-type galaxies that end up with high angular mo¬ 
mentum at z = 0 (top-left panel, selected to have the highest 
j*/M^ 3 values) evolve gradually, approximately along the re¬ 
lation, similarly to late-type galaxies (which are not shown). 
The evolution of early-type galaxies with low angular mo¬ 
mentum at z = 0 (top-right panel, with the lowest j*/M~J 3 val¬ 
ues) is less smooth, and generally has lower, sometimes neg¬ 
ative, slopes on the j*— M* plane. The bottom panels show 
two typical kinds of tracks among those galaxies shown in 
the top-right panel (the circles at the ends of those tracks are 
marked with a special color to distinguish them in the top- 
right panel). The bottom-left presents galaxies that initially 
evolve similarly to galaxies with high angular momentum, but 
at some point experience a sharp specific angular momentum 
drop, by > 0.5 dex ("growth-then-drop"). In two of the three 
featured cases, the drop occurs during a major merger. The 
bottom-right panel presents galaxies that experience an early 
(z > 2) rapid specific angular momentum build-up, with peak 
values reaching the relation for late-type galaxies. Later they 
grow only in mass, by ~ 1 dex, but not in specific angular mo¬ 
mentum, which brings them to low / m]/ 3 values ("arrested 
growth"). These galaxies may or may not experience a further 
late-time merger-associated drop. 

A picture emerges in which all early-type galaxies have 
high angular momentum typical for late-type galaxies at some 
point in their past history. Early-type galaxies form a parallel 
relation to late-type galaxies because they start at early times 
on a similar relation, but later some of them experience pro¬ 
cesses that change their j*/M:J 3 values in a roughly mass- 
independent way. Major mergers are responsible for some of 
the "growth-then-drop" type of evolution. The nature of the 
other prominent evolution type, which features a plateau in 
specific angular momentum, and its possible relation to AGN 


feedback, will be explored in future work. We also note that 
the specific angular momenta of ga laxie s exhibit correlations 
with the spins of their dark halos. iTeklu et all (120151) found 
such correlations in the Magneticum Pathfinder hydrodynam- 
ical cosmological simulations, and they are also present in the 
Illustris simulation. While the focus of the current work is on 
the effects of feedback, future work will explore the relation¬ 
ships between the effects of feedback, formation history, and 
dark matter halo spin on galactic angular momentum. 

4. SUMMARY 

In this paper, we study the angular momentum content of 
galaxies in the Illustris simulation. The combination of vol¬ 
ume and resolution allows us to separate the galaxy popu¬ 
lation into categories based on star-formation rate and mor¬ 
phology, and to compare these to analogous observations and 
to theory predictions. The combination of accurate hydrody¬ 
namics and two channels of feedback, namely galactic winds 
and AGN, results in a diverse galaxy population with approxi¬ 
mately realistic angular momentum content. Our main results 
are as follows. 

• When we split the simulated galaxy population at z = 
0 into two groups based on either their specific star- 
formation rate, flatness, or concentration, we obtain two 
approximately parallel relations between specific angu¬ 
lar momentum and stellar mass, in agreement with ob¬ 
servations. 

• When the angular momentum is considered as a func¬ 
tion of the host halo mass, it is found that simulated 
late-type galaxies have an ‘angular momentum reten¬ 
tion factor’ of rjj ~ 100%, while simulated early-type 
galaxies have only rjj ~ 30% of their original specific 
angular momentum. 

• Early-type galaxies with low angular momentum at 
z = 0 are found at high redshift to lie on the high angular 
momentum relation typical for late-type galaxies. Their 
subsequent evolution in specific angular momentum 
and stellar mass is approximately mass-independent, 
hence the parallel, lower relation between these quanti¬ 
ties on which they lie at z = 0. 

• We also explore a smaller cosmological volume run 
with varying baryonic sub-grid models. We find a posi¬ 
tive correlation between the wind mass-loading factor 
and the resulting angular momentum retention factor 
rjj. Galaxies in a simulation without galactic winds 
have the lowest mean angular momentum, correspond¬ 
ing to rjj ~ 30%. 

• In contrast, AGN-driven thermal ‘radio’ bubbles in our 
simulations generally reduce the specific angular mo¬ 
mentum of galaxies by ~ 20% - 50%. 

Further study is required to determine the exact mechanism 
by which galactic winds negate angular momentum losses and 
drive late-type galaxies to ‘equipartition’ of specific angular 
momentum with their dark matter halos. Future studies will 
also need to explore evolution tracks where early-type galax¬ 
ies lose specific angular momentum, and determine their re¬ 
lation to formation history and AGN feedback. Nevertheless 
we can already see from the results presented here that the 
dominant type of feedback during the formation of a galaxy is 
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Figure 4. Evolution tracks of early-type galaxies in mass and specific angular momentum. The simulated z = 0 galaxies visually-selected (Vogelsberger et all 
12014bf) as late-type (black) and early-type (magenta) appear in all panels. Each panel also shows the evolution tracks of a certain sub-population of early-type 
galaxies, where the curve color indicates redshift. Each track ends at z = 0 with a circle. 


an important driver in setting its angular momentum content, 
with late-type formation relying on galactic wind feedback, 
and early-type formation on AGN feedback. Furthermore, 
feedback introduces correlations between the angular momen¬ 
tum, morphology, concentration, and star-formation activity 
in our simulated galaxies in a manner very reminiscent of the 
Hubble sequence of real galaxies. 
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